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A STUDY OF PLANETARY METEOROLOGY 


by George Ohring 

GCA Corporation, GCA Technology Division, Bedford, Mass. 


SUMMARY 

Two essentially independent pieces of research were completed dur 
ing the contract period. Titles and abstracts of these two research 
efforts follow. 

The Temperature and Ammonia Profile in the Jovian Atmo- 
sphere from Inversion of the Jovian Emission Spectrum 


Inversion techniques developed for remote soundings of the earth's 
atmosphere by meteorological satellites are applied to the few earth based 
observations of the Jovian emission spectrum. With a minimum of assump- 
tions, it is possible to obtain directly from the observed spectra the 
profiles of atmospheric temperature and ammonia abundance. The tempera- 
ture profile is obtained from inversion of the radiances observed in the 
infrared bands of methane and hydrogen. With the use of the derived tem- 
perature profile, the ammonia abundance profile is obtained from radiances 
observed in the infrared and microwave absorption bands of ammonia. The 
temperature profile is characterized by a definite tropopause region with 
a temperature of about 115K and a stratospheric region in which the tempera- 
ture slowly increases with altitude. The derived ammonia profile indicates 
the presence of a saturated ammonia layer with a base temperature of a* 140K. 
It is concluded that inversion of infrared and microwave emission spectra 
provides a powerful technique for sounding the Jovian atmosphere, and that 
more accurate measurements should be obtained from fly-bys and earth based 
observations . 

A Technique to Deduce the Martian Atmospheric Temperature 
and Water Vapor Profiles from Orbiter Observations of 
the Planet's .IR Limb Radiance Profile 

The concept is described of deducing the temperature and constituent 
profile of a planetary atmosphere from orbiter measurements of the planet's 
: IR limb radiance profile. Expressions are derived for the weighting func- 
tions associated with the limb radiance profile for two infrared transmis- 
sion models. Analysis of the weighting functions for the Martian atmosphere 
indicates that a limb radiance profile in the 15p. CO 2 band can be used to 
determine the Martian atmospheric temperature profile from 20 to 60 km. 
Simulation of the Martian limb radiance profile in the rotational water 
vapor band indicates that Martian water vapor mixing ratios can be in- 
ferred from limb radiance observations in a water vapor band. 
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SECTION I 


INTRODUCTION 


The primary objectives of the research contract were: 

(1) To study the meteorology of the planet Jupiter, with particular 
emphasis on the deduction of the atmospheric temperature profile through 
analysis of available infrared and microwave observations of the planet's 
emission spectrum. 

(2) To develop inversion techniques for inferring atmospheric tem- 
perature and constituent profiles of planetary atmospheres. 

Work was performed on both of these topics, although the emphasis was 
placed on item (1) . As a result of the work on item (1) , temperature pro- 
files have been derived for the Jovian atmosphere. In addition to the 
temperature profiles, it was also possible to derive profiles of the 
ammonia abundance. This work is summarized in Section II. The effort on 
item (2) was devoted to development of the equations for inferring the 
temperature and water vapor profiles of planetary atmospheres from orbiter 
observations of the planet's limb radiance profile in C0 2 and H 2 0 absorp- 
tion bands. The feasibility of applying the technique to the Martian 
atmosphere has also been examined. This work is reviewed in Section III. 

The work on both of these topics makes use of inversion techniques 
that are being used or will be used for deducing temperature and constituent 
profiles of the earth s atmosphere from meteorological satellite observa- 
tions. The results indicate that these techniques are very useful for the 
general problem of remote sounding of other planetary atmospheres. 

Efforts should be made to obtain more accurate observations of the 
Jovian emission spectrum in the infrared and microwave, and the feasibility 
of deducing the temperature and water vapor profiles of the Martian and 
Venusian atmospheres should be further examined. It is also suggested 
that the technique that was applied to determine the NHo abundance profile 
on Jupiter, be applied to the determination of the water vapor profile on 
Venus from existing microwave emission observations. 
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SECTION II 


THE TEMPERATURE AND AMMONIA PROFILE IN THE JOVIAN ATMOSPHERE FROM 
INVERSION OF THE JOVIAN EMISSION SPECTRUM 


A. Introduction 

Taylor (1972) has recently suggested an experiment to measure the 
vertical temperature profile of the Jovian atmosphere by fly-by observa- 
tions of the planet* s infrared emission spectrum in the 7,5\x methane band. 
The temperature profile would be deduced from the spectrum with the use 
of inversion techniques that have been developed for satellite tempera- 
ture soundings of the earth* s atmosphere (Chatline, 1970; Smith, 1970). 

There is already available from earth based observations a measured spec- 
trum for Jupiter for the region 2.8 to 14^, which includes the 7.5|j, methane 
band (Gillett et al., 1969). The observed spectrum also includes absorp- 
tion bands due to hydrogen and ammonia. Both hydrogen and methane can be 
assumed to be uniformly mixed in the Jovian atmosphere, and relative 
abundance estimates are available for these gases (see, for example, 

McElroy, 1969). Hence, the observed radiances in both the methane and 
the hydrogen bands provide information on the temperature profile. 

Ammonia is not uniformly mixed in the Jovian atmosphere; it appears 
to be saturated over at least a part of the atmosphere that has been ob- 
served spectroscopically (McElroy, 1969) and, hence, its mixing ratio 
probably varies considerably with height. Thus, unless the ammonia pro- 
file is known a priori, which it is not, observations in the ammonia band 
cannot be used for temperature inferences. But , if one can derive the 
temperature profile from the radiances observed in the methane and hydro- 
gen absorption bands, then, in principle, one can use the radiances observed 
in the ammonia bands to derive an estimate of the ammonia profile. The 
problem is similar to that of deriving the water vapor profile in the 
earth *s atmosphere from satellite observations of an infrared water vapor 
absorption band, given the temperature profile (see, for example, Smith, 
1970). In addition to its absorption bands in the infrared region, ammonia 
has an absorption band in the microwave region centered near 1.25 cm. Earth 
based observation of the Jovian emission in this microwave band are also 
currently available (Law and S tael in, 1968; Wrixon et al., 1971) and can 
be used together with the infrared ammonia bands to estimate the Jovian 
ammonia profile. 

In this report, we attempt to derive estimates of the temperature and 
ammonia profiles in the Jovian atmosphere from earth based observations 
of the planet *s infrared and microwave emission, using inversion tech- 
niques developed for analyzing meteorological satellite spectra of the 
earth *s atmosphere. Admittedly, the quality of the Jovian spectra is not 
as high as the satellite spectra of the earth *s atmosphere. On the other 
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hand, our knowledge of the Jovian atmosphere is not as great as our 
knowledge of the earth's atmosphere, so that the accuracy demands of 
soundings of the Jovian atmosphere are not as severe. 


B. Background 

Much of our information on the composition and structure of the Jov- 
ian atmosphere derives from spectroscopic observations of the planet's 
reflection and emission spectra and theoretical model building. The ob- 
served reflection spectra can be the result of cloud reflection, and 
atmospheric absorption and scattering processes, and are difficult to 
interpret (see, for example, McElroy, 1969). However, it has been pos- 
sible to obtain from such observations estimates of relative and absolute 
abundances of atmospheric gases and rotational temperatures using a re- 
flecting layer analysis. The abundances derived from such an analysis 
refer to the abundances above some reflecting level; the rotational tem- 
peratures refer to some sort of effective temperature of the atmosphere 
above the reflecting level. From observations of the planet's emission 
at infrared and microwave wavelengths additional estimates of atmospheric 
temperatures have been derived. However, the altitudes to which these 
temperatures refer are uncertain. 

Theoretical model building makes use of available observations and 
applicable theory to develop models of the atmosphere consistent with the 
observations. This type of activity has taken several forms. For example, 
Kuiper (1952) and, most recently, Owen (1969) have pieced together are 
atmospheric composition and structure model from observations made at many 
wavelengths together with certain assumptions on the temperature and ammonia 
profiles. Trafton (1967) and Trafton and Munch (1969) have computed radia- 
tive-convective equilibrium temperature profiles based upon observational 
evidence of atmospheric composition. Hogan et al. (1969) have computed 
radiative-convective equilibrium temperature profiles for a number of dif- 
ferent boundary conditions and compared the emission spectrums for these 
profiles with the observed infrared and microwave spectrums. Wildey and 
Trafton (1971) have attempted to match infrared limb darkening observations 
of Jupiter in the 8-14p. region with theoretical limb darkening curves com- 
puted for various assumptions on atmospheric composition and temperature. 
Wrixon et al. (1971) have compared their observed microwave spectrum with 
one computed for an assumed Jovian model atmosphere. 

On the basis of the observational evidence and the results of the 
theoretical model building, the following general picture emerges. The 
Jovian atmosphere consists mainly of hydrogen and helium, with the hydro- 
gen percentage being at least 67 percent, and the hydrogen abundance above 
the reflecting level being about 67 km-atm (McElroy, 1969). Estimates of 
the CH^ abundance range from 30m atm to 150m atm, with a recommended value 
of 45m atm (McElroy, 1969). Ammonia is apparently saturated over part of 
the Jovian atmosphere. Ammonia abundances derived from analyses of the 


4 



2 3 

6470$ band are 7 x 10 cm atm (Kuiper, 1952) and 1.2 x 10 cm atm (Owen, 

1969). Abundances derived from analyses of ultraviolet observations 
yield an upper limit of ~ 2 x 10“^ cm atm, which refers to an effective 
reflecting level higher in the atmosphere (due to Rayleigh scattering) 
than does the abundance derived from the 6740$ band (Anderson et al., 1969) 
There are apparently two cloud layers present in the observable Jovian 
atmosphere (Danielson and Tomasko, 1969; Owen, 1969). The lower cloud ap- 
pears to be dense and at a temperature of ~ 220K with composition unknown; 
the upper cloud appears to be optically thin in the near and far infrared, 
composed of ammonia crystals, and at a temperature of ~ 150K. Pressures 
at the lower cloud-top appear to be of the order of several atmospheres; 
at the upper cloud of the order of a fraction of an atmosphere. 

It is interesting to note that, in much of the theoretical model build 
ing, and also in the interpretation of some of the spectral observations, 
the approach is one of determining which one of a number of different model 
atmospheres best satisfies a particular observation, or set of observations 
This is a rather indirect way of arriving at estimates of the atmospheric 
models most consistent with a given observation. For the analysis of the 
infrared and microwave emission spectra, a more direct technique — inver- 
sion of the actual observations — can be applied to deduce the temperature 
and ammonia profiles as a function of pressure. From the temperature- 
pressure profile, the temperature and ammonia profiles versus relative 
height can be obtained. The techniques for accomplishing such inversions 
are described in the next section. 


C. Inversion Techniques 


Temperature inversion .- The observed radiance at wavenumber for 
an optically thick atmosphere, such as the Jovian atmosphere in the infra- 
red, can be written as 


N^) 


1 



B[v i ,T(p) ]dt (v^p) 


o 


( 1 ) 


where t(v^,p) is the transmittance from the pressure level p to the effec- 
tive top of the atmosphere, and B is the Planck function, given by 


B (v ,T) 


c l v 


3 /[exp(c 2 v/T)-l] 


( 2 ) 


where T is the temperature and c-^ and c 2 are known constants. Our objec- 
tive is to derive the temperature-pressure profile from observations of 
the radiances at n wavenumbers, given the distribution of the absorbing 
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gases, which determines the transmittances . Smith (1970) has derived an 
iterative technique for deriving the temperature profile from such obser- 
vations. One starts with a guess temperature profile and iterates until 
some convergence criterion is achieved. 

In an iterative form, equation (1) becomes 


N(v i ) - 


W 


<v 


/{B j+1 [v i ,T(p)] - B 3 [v i ,T(p)]}dT(v i ,p) 


(3) 


where j and j+1 represent the j th and j+1 iterations. N( v ± ) is the meas- 
ured radiance and Nj(v^) is the radiance computed from equation (1) using 
the temperature profile at step j. For purposes of deriving the iteration 
technique it is assumed that for each wavenumber the difference of Planck 
functions within the integral is independent of pressure within the sensed 
atmospheric layer. Hence, the following iterative equation results 


j+1 -i 

B (v^p) = B J (v i ,p) + [N(v i ) 


N- 


(v i )] 


(4) 


With the use of the Planck formula, it follows that an independent 
estimate of the temperature at the pressure level p can be obtained for 
each observing wavenumber from 


j + 1 / 

T (v^p) = c^/ln 


3 j +1 

C i V i + B (v i'P> 


L B 3+1 (v i# p) 


(5) 


The best new estimate of the temperature at the pressure level p is 
given by the weighted average of the estimates from the n observed wave- 
numbers 


T 


j+1 


(P) 


i+1 11 

= I T J (v i ,p)W(v i ,p)/ Z W(v, ,p) 

i=l i=l 1 


( 6 ) 
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where, from equation (3), it follows that the weights should be 


W(v i# p) = dx(v i# p) 


( 7 ) 


The convergence criterion generally used in analyzing satellite ob- 
servations of the earth's atmosphere is that 


N j+1 (v i ) - N j (v i ) 


/N(v^) <_ e i = l,2,...n 


( 8 ) 


where € is a fraction of a percent. This value of e is commensurate with 
the accuracies of the observations. For our analysis of the Jovian atmos- 
phere spectrum, we used a different convergence criterion. The iterations 
were stopped when changes in temperature from one iteration step to the 
next were < 0.5K at all atmospheric levels. The radiance observations in 
the infrared methane and hydrogen absorption bands are used in the inver- 
sion procedure. 

Ammonia inversion .- We now turn to the problem of deriving the ammonia 
profile from the observed infrared and microwave radiances in the ammonia 
absorption bands. The temperature profile is assumed to be known from the 
temperature inversion procedure. By integrating equation (1) by parts, one 
can obtain 


B[v i# T(p r ) 

N(V ) = B[V ,T(p )] +/ T(Vi,p)dB((v if T(p)] 

B[v i ,T(p 0 ) 


(9) 


where p is the effective top of the atmosphere for ammonia emission and 
p is a reference pressure deep in the atmosphere. In a manner similar to 
that used for deriving the temperature iteration equations, Smith (1970) 
shows that the abundance, a, of optically active gas above a pressure level 
p can be estimated at iteration step j+1 from 


a^ +1 (v i# p) = a3(p){l + [NfVi) - N? (v^ J/S^} 


( 10 ) 
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where 




B[Vi,T(p r )] 

/ (p) 

B [v^ ,T (p Q ) ] 


3t (V i ,p) 


3a(p) 


dB[Vi fT (p)l 


(ID 


The best estimate of the abundance above the pressure level p is 
given by a weighted average of the abundances obtained from (10) 


a^ + 1 (p) 


m 


m 


B E a^ +1 (v.,p) (v,p)/ l W^(v.,p) 
i=l i=l 


( 12 ) 


where the weights are given by 


W (v i# p) 


3x (v i ,p) 
3a (p) 


j 

dB(v i ,p) 


(13) 


and m is the number of wavenumbers for which radiances are observed. 

Equation (13) is appropriate for an absorption band in a single spec- 
tral interval. For observations in which absorption bands in different 
parts of the emission spectrum are used in the inversion, the Planck 
function will vary considerably due to its wavenumber dependence. Hence, 
equation (13) would tend to weight more heavily observations at those wave- 
numbers near the peak of the Planck function at the expense of observations 
at the other wavenumbers. This is not desirable and we have removed this 
distortion by assuming that the weights are given by 

pTCv-.pfP 

- piftoj < 14 > 


The iteration procedure is similar to that for the temperature pro- 
file. One starts with a guess distribution of ammonia abundance and 
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iterates until some convergence criterion is met. As opposed to the 
temperature profile inversion, however, the weights change from iteration 
to iteration since they are themselves determined by the ammonia profile. 
The convergence criterion adopted for the analysis of the Jovian ammonia 
radiances is 


[^(Vi) - (v 4 ) 1/N? (vp 


< 

_ e 


0.01 i » 1,2, . .m 


(15) 


Transmittance models . - To apply the inversion techniques, one needs 
transmittance functions for the absorption bands. For the 7 .7\x band of 
CH^, we have followed McClatchey et al. (1970). McClatchey et al. present 
an empirical transmittance function in graphical form, from which we have 
derived an analytical representation of the transmittance 


x(v,p) = exp {-[k v L(p) ]°* 56 } 


(16) 


where t(v,p) is the transmittance of a 20 cm interval centered at wave- 
number v from the pressure level p to the top of the atmosphere; k v is 
the absorption coefficient for the interval; and L is a pressure corrected 
path length from the level p to the top of the atmosphere 


L> (p) 


1.75 An 


1.75 

P 


(17) 


where n, is the CH, number density, n is the atmospheric number density, 
and 


* / ,0.75 

A = n (p) ug 
STP s 


(18) 


where I s the number density at STP, p g is standard pressure, \x is the 

mean molecular mass of the atmosphere, and g is the gravitational accelera- 
tion. Values of the absorption coefficients derived from this analysis 
are shown in Table 1. Further details may be found in Ohring (1971) • 

For the pressure induced absorption of hydrogen, we follow Trafton 
(1967). The transmittance is given by 
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TABLE 1 


ABSORPTION COEFFICIENTS FOR THE 7.7y BAND OF METHANE 


Wavenumber Wavelength 

(cm - l) (y) (atm-cm)"l 


1160 

8.6 

5.4 

X 

10 -3 

1200 

8.3 

6.2 

X 

10 -3 

1220 

8.2 

1.7 

X 

10~ 2 

1240 

8.1 

1.1 

X 

10" 1 

1260 

7.9 

6.2 

X 

10" 1 

1280 

7.8 

9.4 

X 

10 -1 


10 






( 19 ) 


00 

t(v,») = exp[ -/ k dz] 
z v 


where t(v,z) is the transmittance from the height z to the top of the 
atmosphere, and k y is the volume absorption coefficient of H 2 > which, for 
a hydrogen helium mixture, can be written as 


n 2 

[VT> + qB v (T7f 


( 20 ) 


where n 2 is the hydrogen number density, c is the speed of light, q is 
the ratio of helium to hydrogen number densities, and A V (T) and B V (T) are 
the pressure induced binary absorption coefficients corresponding to ab- 
sorption in pure H 2 and the enhancement in He-H 2 mixtures, respectively, 
at temperature T. Trafton and Munch (1969) and Trafton (1967) present 
graphs of A v and B v as a function of wavenumber for two different tempera- 
tures: T = 160K and T = 80K. From these graphs we have derived the values 

of A and B for use in equation (20). 
v v 

Two ammonia bands are present in the infrared spectrum observed by 
Gillett et al. (1969). For the 10.5fi NH 3 band, we follow Gille and Lee 
(1969), who calculated the positions, half-widths, and intensities of ab- 
sorption lines due to NH 3 at v < 1400 cm" 1 -, and used this information to 
calculate the parameters of the Goody random absorption model with exponen- 
tial distribution of intensities. The transmittance for a spectral inter- 
val, Av, 25 cm "1 wide, and centered at v, from the height z to the top of 
the atmosphere can be written as 



( 21 ) 


where a is the ammonia abundance above level z, P e is the mean effective 
pressure and f the mean effective temperature along the path (for correct- 
ing for the pressure and temperature dependences of the half widths of the 
absorption lines), and values of cr v and P v are the parameters of the Goody 
random model. Values of cr v and p y are tabulated by Gille and Lee (1969) 
for a number of different temperatures. We have neglected the temperature 
dependence, which is relatively small for the temperature and wavenumber 
range of concern in the Jovian atmosphere, and have adopted Gille and Lee*s 
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values of <r v and P v for a temperature of 150K. We have also adopted Gille 
and Lee's expressions for the effective pressure and temperature. 

For the 6.1(i NH3 band we have fit a strong line transmittance law to 
the laboratory observations of absorption by France and Williams (1966). 
The transmittance in the strong line approximation can be written as 


_ _ 1/2 

x v - exp { - k^a P g f] ) 


( 22 ) 


It is assumed that P g and f have the same forms as for the 10. 5p NH, band. 
Values of k y derived from this analysis are tabulated in Table 2. Unfor- 
tunately, the values of k derived from this analysis must be considered 
very uncertain since the laboratory data include only one set of observa- 
tions for NH3 paths as large as 49 atm-cm, the remaining observations being 
for paths less than 3 atm-cm. 

For absorption by ammonia in the microwave region, we follow Wrixon 
et al. (1971). In general, the ammonia absorption coefficient is given by 
(Townes and Schawlow, 1955) 


k 


v 


00 

z 

J=1 


J 

z 

K=1 


1.225P 



x 


S ( K) I2J+1)K 2 
J(J+1) 


exp { - [2 .98J (J+l) - 1 . 09 K 2 ] 4 .8/T}F (v) 


(23) 


where k^ is the absorption coefficient at frequency v in units of cm” 

(note that for the microwave region v represents frequency, not wavenumber), 
Pq is the partial pressure of NH3 in atmospheres, S(K) = 3 for K a multiple 
of 3, and S (K) = 1.5 otherwise, and F(v) is the line shape factor. At low 
pressures (p < 10 cm Hg of pure NH3) , the absorption lines are distinct) 
and individual lines can be adequately represented by the Van Vleck-Weisskopf 
shape. As the pressure is increased, the individual lines are broadened and 
wiped out, and a single Van Vleck-Weisskopf expression can be used to express 
the absorption coefficient. As the pressure is further increased (p > 30 cm 
Hg of pure NH^), the shape of the absorption spectrum changes from the res- 
onant Van Vleck-Weisskopf form to the non-resonant Debye line shape (Bleaney 
and Loubser, 1950). Ben Reuven (1966) has derived a formula that approxi- 
mates the absorption at these pressures. In the present study, we use the 
Van Vleck-Weisskopf shape — summing the individual lines as indicated by 
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TABLE 2 


ABSORPTION 

COEFFICIENTS FOR THE 

6. In BAND OF AMMONIA 

Wavenumber 

(cm’l) 

Wavelength 

(n) 

(cm- atm z ) 

1785 

5.6 

0.10 

1850 

5.4 

0.023 


equation (23) — at the upper levels (lower pressures) of the Jovian atmo- 
sphere, and Ben Reuven's formulation at the higher pressures. The change- 
over is made when the line widths in the Jovian atmosphere are equal to or 
greater than the line widths at a pressure of pure NH 3 of 30 cm Hg. The 
expressions for the Van Vleck-Weisskopf and Ben Reuven line shape factors 
given by Wrixon et al. (1971) were adopted for this study. 

The transmittance at frequency v from the level z to the top of the 
atmosphere is then given by 



(24) 


D. The Observations 

The infrared spectrum used in the present study is based upon the 
Gillett et al. (1969) observations of Jupiter over a period of five dif- 
ferent nights. Observations were made in the 2.8-5.6p, and the 7.4-15p 
regions with a resolution of AX/X 0.02. The observations were combined 
into the single spectrum shown in Figure 1. Statistical errors in the ob- 
servations can be judged by the vertical spread of adjacent points. 

For X < 4|i, Gillett et al. assume that reflected and scattered sun- 
light is the primary source of the radiation; for X > 5p they assume that 
thermal emission from Jupiter is the primary source. The shape of the 
spectrum at X > 5p can be attributed to ammonia bands in the regions 
5.2 - 5 . 6|_i and 9 - 12 p; to a methane band in the region 7.8 - 8.4p; and 
to hydrogen absorption at X > 12p, (Gillett et al. 1969). We have smoothed 
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Figure 1. Observed Jovian infrared emission spectrum (after 
Gillett et al., 1969). 


the spectrum by eye and read off the emission temperature at a number of 
wavelengths in each of these spectral regions. The uncertainty in reading 
off the temperatures may be as large as 2K. The observed emission tempera- 
tures are shown in Table 3. Observations in the CH^ and H 2 bands are used 
to determine the temperature profile; observations in the NH 3 bands are 
used to determine the ammonia profile. 

The emission spectrum of Jupiter in the vicinity of the 1.25-cm 
ammonia band at microwave wavelengths has been observed by Law and S tael in 
(1968) and Wrixon et al. (1971)1 The measurements of Law and Staelin cover 
five wavelengths between 25.4 and 19.0 GHz; those of Wrixon et al. cover 
eight wavelengths between 20.5 and 35.5 GHz. Additional observations at 
individual wavelengths within this region have been obtained by Kellerman 
(1970) at 15.8 GHz, Tolbert (1966) and Kalaghan and Wulfsberg (1968) at 
34.9 GHz and Thornton and Welch (1963) at 36 GHz. The measured brightness 
temperature and their associated error flags are plotted in Figure 2, after 
Wrixon et al. (1971). As can be seen, there is quite a bit of uncertainty 
in the shape of the spectrum. We have attempted to draw a smooth curve by 
eye through the data on the high frequency side of the absorption band, 
giving more weight to the Wrixon et al. brightness temperatures because of 
their smaller error bars. Temperatures at three frequencies were then read 
off, to be used, along with the infrared NH 3 band observations, in the in- 
version for the ammonia profile. These microwave temperatures are shown 
in Table 4. 


E. Results 


Temperature . - The entire concept of remote probing by measuring the 
emission spectrum is based upon the fact that the observation radiation 
originates in different parts of the atmosphere, depending upon the ab- 
sorption coefficient, and hence, the wavelength of the radiation. This 
can be seen from Figure 3, which is a plot of weighting functions, dx/dz, 
for the CH^ and H 2 wavelengths used in this study. The weighting func- 
tions are proportional to the weights of equation (7) and indicate where 
in the atmosphere the radiation at each wavelength originates. These 
weighting functions are for an atmosphere consisting of the following 
volumetric fractions: H 2 = 0.67, He = 0.33, CH^ = 4.7 x 10"^. As in all 

the computations discussed in this paper, a vertical spacing of 2 km is 
used in the numerical integrations, and the zero reference level for alti- 
tude is taken to be the pressure level 5 atm. 


Figure 3 indicates that at these wavelengths and for this composition, 
information is received on the temperature profile in the height range 
25 to 100 km (if we assume that no information is received for altitudes 
below the altitude of the half-peak value of the lowest weighting function, 
or for altitudes above the altitude of the half-peak value of the highest 
weighting function). The hydrogen band weighting functions are narrower 
than the methane weighting functions due to the second power dependence of 
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. TABLE 3 


OBSERVED JOVIAN EMISSION TEMPERATURES (T ) FOR 

£ 

THE INFRARED BANDS OF CH 4 , H 2 , AND NH 3 
(From Figure 1) 


Absorption Band 

Wavenumber 

(cm -1 ) 

Wavelength 

(y) 

t e 

(K) 

CH 4 

1160 

8.6 

138 


1200 

8.3 

138 


1220 

8.2 

138 


1240 

8.1 

139 


1260 

7.9 

141 


1280 

7.8 

145 

H 2 

715 

14.0 

120 


800 

12.5 

125 

NH 3 

862 

11.6 

127 


937 

10.7 

125 


1037 

9.6 

125 


1785 

5.6 

170 


1850 

5.4 

190 
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JUPITER 



Figure 2. Observed Jovian microwave emission spectrum in the region of the 1.25-cm ammonia band. 

Error flags with solid lines are observations of Wrixon et al. (1971); error flags with 
dashed lines are observations of Law and Staelin (1968) and others. Curve through obser- 
vations on high frequency side represents smoothing of the observations to derive bright- 
ness temperatures for use in the inversion procedure. 



TABLE 4 


OBSERVED JOVIAN BRIGHTNESS TEMPERATURES, T 

B' 

FOR THE NH 3 MICROWAVE BAND 
(From Figure 2) 


Frequency 

Wavelength 

T 

(GHz) 

(cm) 

(K? 

35 

0.86 

145 

30 

1.0 

135 

24 

1.2 

130 
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Height (Km) 


Wavenumber » (cwrj) 



Figure 3. Weighting functions for the infrared absorp 

tion bands of CH, and H 0 . 

4 2 


0.06 
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the hydrogen absorption coefficient on the hydrogen density. It is of 
interest to note that, for this atmospheric model, radiances in the two 
hydrogen wavenumbers and in the 1160 and 1200 cm“l methane wavenumbers 
provide information on the atmospheric temperature in approximately the 
same altitude region. 


We have computed the temperature profiles for several atmospheric 
composition models, covering the range of present uncertainties in the 
hydrogen and methane amounts, as discussed in Section B. The hydrogen, 
helium, and methane fractional volumes for these models are summarized in 
Table 5. In all calculations, the temperature profile was extended adia- 
batically from the lowest point at which information is received from the 
temperature inversion to the reference pressure of 5 atm. 

Temperature profiles for the first set of experiments, in which the 
H2 fractional volume is assumed to be 0.67, are shown as curves (a) and (b) 
in Figure 4. The temperature profile for a CH4 fractional volume of 
3.1 x 10“^ is almost identical with that for a CH4 fractional volume of 
4.7 x 10“^, and, hence, is not shown. The differences in the temperature 
profiles are due to a slight shift upwards of the CH4 weighting functions 
with respect to the H2 weighting functions with the higher CH4 amount. 
Basically, however, the profiles are similar, indicating that the derived 
temperature profile is not particularly sensitive to the present uncertainty 
in CH^ amount. The temperature profile is characterized by a definite tro- 
popause and stratospheric region similar to that in the earth's atmosphere. 
Gillett et al. (1969), interpreting their observed spectrum on the basis of 
a semiquant itative analysis, deduced the shape of the temperature profile 
shown in Figure 4 and suggested that solar heating in the 3020 cm"l CH^ 
band is the cause of the temperature increase above the tropopause. Radia- 
tive equilibrium temperature calculations, including absorption of solar 
energy in this band, have reproduced this structure (Hogan et al., 1969). 
Thus, it appears reasonable to conclude that the stratospheric temperature 
structure shown in Figure 4 is indeed real and is due to CH4 absorption of 
solar radiation. 


The RMS percentage difference between radiances computed from the 
derived temperature profiles and the observed radiances is 16 to 17 per- 
cent. One of the reasons for these differences is the relatively large 
differences in emission temperatures observed in the H2 band and those 
observed in the 1160 to 1200 cm -1 region of the CH4 band. As indicated 
earlier, the radiances at these wavenumbers refer to approximately the 
same atmospheric region, and hence, they should yield the same emission 
temperature. The fact that they do not suggests that they are inconsistent. 
The iterative inversion scheme handles this inconsistency by objectively 
weighting each of the radiances so that the resulting temperature is based 
upon a suitable average of the observations. It should also be pointed out 
that in the region of the methane band at 8|i, and at the temperatures of 
the Jovian atmosphere, a one degree change in atmospheric temperature is 
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TABLE 5 


ATMOSPHERIC COMPOSITION (FRACTIONAL VOLUME) MODELS USED FOR 
CALCULATIONS OF TEMPERATURE PROFILE 


Experiment No. 

H 2 

He 

CH, 

4 

i 

0.67 

0.33 

4.7 

x 10‘ 4 

2 

0.67 

0.33 

3.1 

x 10' 4 

3 

0.67 

0.33 

1.6 

-3 

x 10 

4 

1.00 

0 

7.0 

-4 

X 10 

5 

1.00 

0 

4.7 

x 10' 4 

6 

1.00 

0 

2.3 

x 10" 3 
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Height (Km) 



Temperature (K) 

Figure 4. Inferred temperature profiles: (a) Experiment 
2, (b) Experiment 3, (c) Experiment 5. Dashed 
lines represent adiabatic extensions of inferred 
temperature profiles to a pressure level of 5 
atmospheres, which serves as the zero reference 
for the altitude scale. 


equivalent to a 10 percent change in radiance. Thus, 20 percent RMS 
differences between observed and computed radiances in the methane band 
represent temperature differences of only 2K. 

Temperature profiles for the second set of model atmospheres (H 2 = 
1 . 0 ) were almost identical and, hence, only one of this set is plotted 
in Figure 4, as curve (c). The main effect of a pure H 2 atmosphere is 
an increase in the scale height which leads to an increase in the alti- 
tudes of and a broadening of the weighting functions. The derived tem- 
perature profiles for these cases are characterized by much higher tropo- 
pauses. Again, there is little difference in temperature profile with 
varying amounts of CH 4 . RMS percentage differences between radiances 
computed with these temperature profiles and the observed radiances are 
in the range 15 to 17 percent. 

Ammonia.- Making use of the derived temperature profiles, we have 
computed ammonia profiles for the same model atmospheres of Table 5. 
Ammonia weighting functions for the initial guess ammonia profile are 
shown in Figure 5. This ammonia profile is characterized by a total abun- 
dance above the 5 atm level of 1.3 x 10 3 atm cm and a scale height of 
1.6 km. The chosen guess scale height is intended to simulate an NH 3 
saturation condition. The hydrogen volume fraction is 0.67. The NH^ 
weighting functions are narrower than the CH 4 or H 2 weighting functions 
due to the small NH^ scale height. 

Ammonia profiles for the first set of experiments (H 2 = 0.67) are all 
similar; hence, only one is plotted in Figure 6 , that for experiment 2. 

The final NH 3 weighting functions for this ammonia profile are shown in 
Figure 7. It is assumed that no information is obtained on the NH^ pro- 
file below the peak of the lowest weighting function or above the peak of 
the highest weighting function. The quantity plotted in Figure 6 is the 
NH 3 abundance as a function of altitude. NH 3 partial pressures, deter- 
mined from the abundance profile, were compared to saturation vapor pres- 
sures for NH~. The region where NH 3 saturation was found is indicated in 
the graph. Thus, at an altitude of about 30 km, an ammonia ice crystal 
cloud would presumably be present with a base temperature of 142K and 
top temperature of 130K. 

For clarity, in the graph of the NH^ weighting functions associated 
with the NH 3 profile, only one weighting function from each of the three 
NH 3 absorption bands is plotted. The multiple peaks of the final Weight- 
ing functions is due to the uhsmooth character of the final NH 3 profile. 
The unsmooth final profile may, in turn, be due to inconsistencies in the 
observed radiances in the three NH 3 absorption bands. 

Ammonia profiles for the second set of experiments (HL = 1.00) are 
also almost identical; hence, only the computed profile for experiment 5 
is shown (Figure 8 ). The shape of the profile is similar to that of the 
H = 0.67 case. However, the ammonia is now at a higher altitude and its 
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Figure 5. NH 3 weighting functions for the initial "guess" NH profile 
terized by a total abundance of 1.3 x 10 3 atm-cm and a seal 




of 0.67 (Experiment 2). 
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Figure 7. NH^ weighting functions for inferred NH 3 abundance 
profile of Figure 6 . 



I0* 2 I0* 1 10° 10' I0 2 I0 3 (0 4 

NHj Abundance (atm-cm) 


Figure 8 . Inferred NH 3 abundance- profile for, an H 2 
fractional volume of 1,0 (Experiment 5). 



scale height is greater. Both of these differences are ultimately due 
to the greater scale heights of the ^ = 1.00 atmosphere. A saturated 
layer is also found for this set of experiments at an altitude of about 
60 km, with a base temperature of 137K and top temperature of 129K. RMS 
percentage differences between the radiances computed with the final NH^ 
profile and the observed radiances are % 6 percent for this set; for the 
first set, the RMS percentage differences are ** 7 percent. 

The derived NH~ abundances may be compared to some of the values 
found in the literature. With the use of a reflecting layer hypothesis, 
total abundances of 7 x 10^ (Kuiper, 1952) and 1.2 x 1C>3 cm atm (Owen, 
1969) have been obtained from the 6470& observations. Our maximum value 
of abundance is of the order of 4 x 10^ atm cm. This would be in rough 
agreement with Owen's value if the reflecting level is close to the level 
of our maximum value of abundance. The model builders suggest NH 3 cloud 
base temperatures of 150K (Danielson and Tomasko, 1971) and 160K (Owen, 
1969). NH^ cloud base temperatures derived from the present analyses' 
are in the vicinity of 140K. The present analysis supports Danielson and 
Tomasko's suggestion of an upper, thin NH^ ice crystal cloud in the Jovian 
atmosphere. 


F. Conclusions 

Inversion techniques for temperature and atmospheric composition, 
originally developed for analyzing meteorological satellite observations 
of the earth's infrared emission spectrum, have been successfully applied 
to earth based observations of the Jovian infrared and microwave emission 
spectrum. Temperature and ammonia profiles for the Jovian atmosphere have 
been derived. The derived Jovain temperature profile is characterized by 
a tropopause region with a temperature of 115K located in the pressure 
range 0.3 to 0.6 atm and a stratospheric region in which the temperature 
increases with altitude. The actual tropopause temperature may be some- 
what lower than 115K since the temperatures derived from the inversion 
represent essentially vertically averaged temperatures, the vertical reso- 
lution being related to the widths of the weighting functions. The de- 
rived ammonia profile is characterized by a total abundance of % 4 x leP 
atm-cm above the pressure level 2 atm, a rapid decrease of NH^ with alti- 
tude (compared to the rate of decrease of total atmospheric density), and 
a relatively thin (x 5 km) saturated layer at a base temperature of ^ 140K. 

There are a number of uncertainties in the present inversions. These 
include the errors in the radiance observations (particularly in the micro- 
wave region) and uncertainties in the transmittance models used (particu- 
larly for the 6.l[x band of NHg). However, the results are consistent with 
present ideas on the Jovian atmosphere and suggest that inversion of the 
Jovian emission spectrum is a powerful method for sounding the Jovian atmo- 
sphere. With a minimum of assumptions, one can derive temperature and 
ammonia profiles directly from the radiance observations. Improved inver- 
sions can be obtained from fly-by or orbiter observations of the Jovian 
infrared and microwave spectrum, or from more accurate earth based obser- 
vations, together with more accurate transmittance models. 
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SECTION III 

A TECHNIQUE TO DEDUCE THE MARTIAN ATMOSPHERIC TEMPERATURE AND WATER 
VAPOR PROFILES FROM ORB ITER OBSERVATIONS OF THE PLANET'S IR 

LIMB RADIANCE PROFILE 


A. Introduction 

The inference of atmospheric temperature and constituent profiles 
from measurements of the spectral distribution of the emitted thermal 
radiation from a planetary atmosphere has received a good deal of atten- 
tion recently. Such experiments have already been successfully flown in 
the Earth orbiting NIMBUS satellites (see, for example, Conrath et al., 
1970), have been flown on the Mariner 9 flight to Mars (Hanel et al., 
1972), and will be flown on future flights to Mars and Venus. Another 
technique, based upon scanning the limb of the planet rather than scan- 
ning in wavelength, can also provide valuable information on the vertical 
profiles of temperature and constituents of a planetary atmosphere. This 
concept is just beginning to receive attention (Gille and House, 1971) 
and will be tested on the NIMBUS F meteorological satellite. This tech- 
nique has certain advantages over the spectral scan, and can be used to 
complement spectral scans . 

The basic geometry of the limb, or horizon, radiance observation is 
shown in Figure 1. A high angular resolution (~ 0.1 degrees for orbital 
altitude of 10^ km) radiometer on a planetary orbiter intercepts radia- 
tion emanating from an atmospheric path that is tangent to the planetary 
surface at height h, which is termed the tangent height. In the figure, 

R represents altitude. The wavelength interval to which the radiometer 
is sensitive is determined by the nature of the experiment. If the temp- 
erature profile is the object of the experiment, the wavelength interval 
will lie in an atmospheric absorption band of a constituent whose distri- 
bution is known - e.g., the 15p. CO 2 band would be a good choice for Earth, 
and possibly, Mars and Venus. If the water vapor profile is the object 
of the measurement , then an infrared water vapor band would be selected. 
The radiation received by the instrument is a function of the distribution 
of temperature and absorbing gas along the ray path. If the distribution 
of absorbing gas is known, the temperature above the tangent height can be 
inferred from the radiance observations. If the temperature distribution 
is known, we can infer the distribution of an absorbing gas. By scanning 
in tangent height from top to base of the atmosphere, it is thus possible 
to infer the vertical temperature structure or vertical distribution of 
absorbing gas. The radiance versus tangent height curve obtained through 
such a scan is called a horizon radiance profile or limb radiance profile. 
The geometry of the horizon radiance observation technique insures that 
one is receiving information on the temperature, for example, from the 
atmosphere above a particular level — the tangent height. It will be 
shown below that for a gas with constant mixing ratio — such as CO 2 in the 
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Figure 1. Limb Viewing Geometry 


atmosphere of Earth and probably in the atmospheres of Mars and Venus 
(at least to great heights) — the weighting functions for the horizero 
radiance observation are sharply peaked at the tangent height. This tteahs 
that not only is the vertical slice of atmosphere under view limited at 
its lower bound by the tangent height being observed, but its upper bound 
is also effectively limited by the shape of the weighting function. In 
other words, for a particular tangent height h, most of the radiation comes 
from a rather narrow vertical interval between h and h + AZ. Hence * the 
mean temperature of this interval could easily be inferred by an ‘instant 
invers ion 1 technique . 

Below we develop the technique for two infrared transmission models: 
the Goody random model, and the transmission models used by Bartkb and 
Hanel (1968) for calculations of temperature in the upper atmosphere of 
Venus. We then compute weighting functions for the atmosphere of Mars 
using the Goody random model. 


B. Theory 

The intensity of radiation received by an infrared radiometer with 
narrow field of view whose line of sight is tangent to the pianetary 
atmosphere at the tangent height, h of Figure 1, is given by 


I 


oo 



-00 


B ^ d * 




dx 


( 1 ) 


where B is the Planck intensity integrated over the spectral interval Of 
the radiometer, and t is the spectral transmittance of the atmosphere 
from the point x to the radiometer. From the geometry, the following 
relationship between x and height z exists 


2 2 
(R + h) + x 


(R + z) 2 


( 2 ) 


The terms h and z can be neglected in comparison with the other 
terms — since we are interested in the lower part of the atmosphere 
(Z < 100 km) — so that 



(3) 
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and 



( 4 ) 


Thus Equation (1) can also be written as 


I 



(5) 


Expressions will be derived first 


for the weighting functions 


dr 

dz 


which determine where in the atmosphere the radiation received by the sen- 
sor originates. We assume that the observations are made in a COo band 
for purposes of temperature inferences. 

, P ood y Random Model .- If it is assumed that the transmittance is given 
by the Goody random model, then 



( 6 ) 


where k is the mean line intensity, 5 is the mean line spacing, m is the 
effective amount of absorbing gas in g cm -2 along the path, and U is the 
effective half-width along the path. For the Martian and Venusian CO, 
amounts and pressures, and the long paths of concern here, km/^a is much 
greater than unity. Thus, Equation (6) can be written as 
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Following Rodgers and Walshaw (1966), the effective absorber amount 
and half-width are 


m 


$(©)dm 


( 8 ) 


and 



(9) 


where is a correction for the effect of temperature on the line in- 

tensities and is given by 


log e $(e) * a(e - 260) + b(e - 260) 2 , 


( 10 ) 


where Q is temperature and a and b are constants for the band. For 
the present purpose of determining the C0 2 transmittances , an isothermal 
atmosphere is assumed so that 


m a = $(e) 


/“ 


dm 


(in 


S ince 




( 12 ) 


where p is 
and pressure, 


pressure and the subscript 
and dm = 

m co 2 


s 


refers to standard temperature 
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*-p .(f) 


1/2 


(z - h) 


1/2 


dz, 


(13) 


where p C0 o iS the carbon dioxide density, p is the atmospheric density 
and c 2 is the C0 2 mass fraction, Equation (11) can be written 


as 


1/2 1/2 

i5 ‘* (#)o •(%)(»*) CO c / 


p p 


JLl 


(* - h) 


1/2 ( 14 > 


For an isothermal atmosphere both pressure and density have the same scale 
height H. Thus, Equation (14) becomes 


m a * 0(0) a 


/•p \ /e \ 1 ^ 2 / v 1 
(-2) f-i) / R ^ 

VP 8 / \ej \2 J 


1/2 

c P, 


-2 z/H dz 

\Tz-h 


(15) 


where the subscript o refers to surface values on the planet. 


With the substitution, 


2 = 2 (* - h) 
H 


the integral in Equation (15) becomes 


-2h/H 

e 


(2H) 1/2 J 



( 16 ) 


Hie integral in Equation (16) is equal to 2 ^ 1 

erf is the error function, and a = «/2(z-h)/H. Thus 


erf a] 

, where 
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where the negative sign in front of the error function is to be used for 
anterior values of z (to the right of the y axis in Figure 1) and the posi- 
tive sign for posterior values of z. The expression for the transmittance 
from the point z to z = oo can now be written as 



( 18 ) 



Differentiating with respect to z yields 



where the upper signs are to be used for anterior values of z and the 
lower signs for posterior values of z. 
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Bartko-Hanel transmission model .- The Bartko-Hanel (1968) transmit- 
tance for a spectral interval can be written as 


T 


exp 


[- (mu*)"] 


( 21 ) 


where m and n are constants for the spectral interval, and u , the effec- 
tive absorber amount, is given by 


* 


u 




( 22 ) 


where u is path- length of absorbing gas in units of atmo-cm of CC> 2 > T g 
and p s are standard temperature and pressure, respectively, and 7 is a 
constant for a spectral interval. Expression (22) is for a constant 
temperature and pressure path. We assume an isothermal atmosphere and, 
hence, a constant temperature for purposes of evaluating transmittances . 

The pressure does vary over the path and in this case expression (22) pk 
is replaced by p^ , where the bar represents an average over the atmospheric 
path . 


k 

P 



(23) 


From the definition of the path length, 


CCL 


du = 


dx = 


< 5 > 


1/2 


dz 


C0„ 


CO. 


(z " h) 


1/2 


(24) 


where the prime represents conditions of standard temperature and pres- 
sure, and c is the CO 2 mass fraction. 
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The effective absorber amount can then be written as 


u = f (T) p. 


I 


k , 

p du 


( 25 ) 


where 


f(T) 



( 26 ) 


or 


where 


* _ P k 

u = B / p p 


(z 


dz 

- h ) 1/2 


( 27 ) 
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f(T) Pg ' k c 


COr 



( 28 ) 


Under the assumption of an isothermal atmosphere (for transmission 
calculations only), both pressure and density have the same scale height 
H. Thus, 


* 

u = B 


P Q k P Q J ex P [ _ ( k + D z/h| - 

V ( 


dz 


z - h) 


172 ( 2 9 ) 
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where the subscript o refers to the planet's surface. 


If we let 


2 _ (k + 1) (z - h) 
H 


( 30 ) 


the integral in (29) becomes 


(dh) 


1/2 


e' t2 dt 


exp -<k + l)h/H J 

A f(z-h) (k+1) 

“ H 


(31) 


The integral in Equation (31) is equal to O " er ^ , where 


a = v (z-h) (k+1) /H 


(32) 


Thus , 


* k 

u = B p p 
o K o 


1/2 exp [-(k+l)h/H] (1+ erf a) ( 33 ) 


where the negative sign in front of the error function is to be used for 
anterior values of z, and the positive sign for posterior values of z. 
The expression for the transmittance from the point z to z = oo along the 
tangent path can be written as 


t = exp 


mE exp [-(k+l)h/H] (1 + erf a) j 


( n 


(34) 
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where 



( 35 ) 


Differentiating with respect to z yields 


dt 

dz 


(z ’ h ' ) = + n [^^y] / ex P - {raE exp- [(k+l)h/ll] (1 + erf a) } 

, ■ r -1 \ n ~ 1 -a 2 (36) 

|m E exp [-(k+l)h/Hj (1 + erf a)j e 


C. Application to the Problem of Inferring the 
Martian Temperature Profile 

The weighting function is defined as 


W(z,h) 



(37) 


with A and P standing for anterior and posterior, respectively. 

We have used the random transmission model expression (Equation 20) 
to compute weighting functions for the 15|a CO 2 band for the Martian atmosphere 
The weighting functions for tangent heights of 0 to 60 km, at 10 km inter- 
vals in tangent height, are plotted in Figure 2. In these computations it 
has been assumed that the Martian scale height is 10 km, the surface pres- 
sure is 5 mb, the surface density is 1.2xl0“5 g cm“^, the temperature of the 
Martian atmosphere is 200° and the band constants are those given by Rodgers 
and Walshaw (1966). It is obvious from Figure 2 that for tangent heights 
greater than 20 km most of the energy emanates from a rather narrow vertical 
thickness of atmosphere, of order of 5 km. Thus, observations of 15(i band 
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radiance at 5 km or greater tangent height intervals are almost indepen- 
dent of each other and could be used to infer, quite directly, mean 
temperatures with 5 km vertical resolution. A rather simple technique, 
"instant inversion," is developed below. 

From Equation (1) or (5) it can be seen that the intensity of radia- 
tion received by the radiometer can be written as 

1 

1(h) = J B <*t (38) 

t(- °°»h) 


where r(-«>,h) represents the transmittance of the entire atmospheric path 
tangent to the Martian atmosphere at height z = h. The weighting functions 
for tangent heights greater than 20 km are such that most of the energy 
emanates from a rather narrow vertical interval immediately above the 
tangent height. Thus, the temperatures above this narrow vertical inter- 
val do not contribute significantly to the observed emission. It may, 
therefore, _be assumed that B in the integral of Equation (38) may be approx- 
imated by B (h) , which represents the integrated value of the Planck func- 
tion for the average temperature of the layer immediately above the tan- 
gent height. Thus, 


1(h) - B(h) [ 1 - t( - 00 , h)] 


(39) 


term [1 - t(-°°, h) ] may be called the emittance, e (h) • Thus, one can 
infer B(h) from 


B (h) 


1(h) 

€(h) 


(40) 


It is a simple matter to reclaim T(h) from the integrated Planck func 
tion. Hence, the vertical distribution of temperature of the Martian 
atmosphere can be obtained from observations of I as a function of tangent 
height h. 

From Equation (18) 


t(-«>, h) = exp (- V2 Ae~ h/ ^ 
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The emittance of the Martian atmosphere as a function of tangent height 
is shown in Figure 3. 

For purposes of deriving and illustrating the method, we have assumed 
a model isothermal atmosphere to calculate the variation of emittance with 
tangent height. Although the emittance varies with temperature, this is 
a second order effect compared to the variation of the Planck function with 
temperature. Hence, Equation (40) can be used to infer temperatures at dif- 
ferent levels even though e(h) is computed for an isothermal atmosphere. 

The effect of uncertainties in the emittance profile upon inferred tempera- 
tures can easily be examined. 

Preliminary analysis of experiment feasibility indicates that it is 
within the state of the art. The desired vertical resolution at the tan- 
gent height is about 5 km. For a Martian orbiter at an altitude of 1000 km, 
this vertical resolution can be achieved with 0.1 degree vertical field of 
view. A similar successful measurement of the Earth*s horizon radiance 
profile in the 15^ CO 2 band (615 cm"! to 715 cm”l) from a suborbital rocket 
flight utilized a radiometer with 0.025 degree vertical field of view 
(McKee et al., 1967). Expected Martian radiances in the 15 \i band are of 
the order of 1/10 of Earth radiances (assuming stratospheric temperatures 
of 150K on Mars and 230K on Earth). The wider field of view suggested for 
the Martian experiment helps to compensate for the lower signal levels. 

Scan times of the order of 1 second for a complete vertical scan, as in 
the successful Earth experiment, are within the state of the art. The 
effective horizontal resolution of an individual temperature inference 
would be of the order of 200 km. 

This analysis shows that the orbiter based horizon radiance technique 
is potentially a powerful tool for measuring the vertical temperature 
structure of the Martian atmosphere above heights of about 20 km. It would 
thus complement an orbiter infrared spectrometer experiment at 4.3p, or at 
15p which could be used to infer temperature profiles for the lowest 20 km 
of the Martian atmosphere. If one wanted to observe the temperature pro- 
file of the lowest 20 km of the Martian atmosphere with the limb radiance 
technique, one has only to slide the spectral interval of the radiometer to 
a region of slightly lower CO 2 absorption. More sophisticated inversion 
techniques can also be used. 


D. Application to the Problem of Inferring Martian 
Water Vapor Mixing Ratios 

For a water vapor absorption band, we may write the following expres- 
sion for the observed radiance, 1(h), at tangent height h 


1(h) = B(h) [1 - T ( -00 > h)] 


(42) 
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Figure 3. Emittance of the 15p. CO 2 band versus tangent 
height for the Martian atmosphere. 
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where B(h) is the integrated value of the Planck function over the water 
vapor spectral interval for the average temperature of the layer immed- 
iately above the tangent height, and h) represents the transmittance 

in the H 2 O band of the entire atmospheric path tangent to the Martian atmo- 
sphere at height z = h. The transmittance may be written as 


t(-“, 


h) <= exp 


(• 


v2 A e 



(43) 


where 



k is the mean line intensity, 6 is the mean line spacing, 0(0) is a cor- 
rection for the effect of temperature on the line intensities, the subscript 
s refers to standard temperature and pressure conditions, the subscript o 
to Martian surface conditions, a is the mean half-width of the absorption 
lines, p is pressure, 0 is temperature, p is density, R is the radius of the 
planet, H is the scale height, and w is the water vapor mixing ratio. This 
equation is based upon the Goody random model, strong line approximation, of 
transmittance (see Rodgers and Walshaw, 1966, whose nomenclature and numerical 
values are used herein). 


We have calculated the radiances, 1(h), for the spectral interval 
160 cm"* to 280 cm"* in the 1^0 rotational band for different water vapor 
mixing ratios, w. For these calculations, it is assumed that the surface 
pressure was 5 mb, the scale height was 10 km, and the atmospheric tempera- 
ture was 200°K. The results are shown in Figure 4. The sensitivity of the 
radiance to the H 2 O mixing ratio is clearly illustrated. 


These results suggest that, if the temperature profile is known (for 
example, from limb radiance observations in the CO 2 band), it would be 
possible to obtain information on Martian water vapor from orbiter obser- 
vations of the limb radiance profile in a water vapor absorption band. 
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